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In supersymmetric theories a field can develop a vacuum expectation value M ^ 10^ GeV, even 
though its mass m is of order 10^ to 10^ GeV. The finite temperature in the early Universe can 
hold such a field at zero, corresponding to a false vacuum with energy density Vo ~ rn^M^ . When 
the temperature falls below Vq^*, the thermal energy density becomes negligible and an era of 
thermal inflation begins. It ends when the field rolls away from zero at a temperature of order 

1/^ ' m, corresponding to of order 10 e-folds of inflation which does not affect the density perturbation 

On . generated during ordinary inflation. Thermal inflation can solve the Polonyi/moduli problem if M 

On ' is within one or two orders of mag nitude of 10^^ GeV. 

^ ■ 

o ■ 

«D ■ I. INTRODUCTION 

Q: 

ly^ I There is at present a 'standard model' of the Universe before nucleosynthesis, which is described in many reviews 
and several textbooks. According to this model, an early era of inflation sets the initial conditions for a Hot Big Bang, 
which starts far above the critical temperature for the electroweak transition (T ~ 100 GeV) and continues without 
^ interruption until the present matter dominated era begins. 

This picture is pleasingly simple, but it is by no means mandatory in the context of current thinking about the 
fundamental interactions beyond the Standard Model. To be precise, it will not be valid if one or more scalar fields 
have a sufficiently large vev (vacuum expectation value) while at the same time having an almost flat potential. 
I The reason is that the particle species corresponding to the oscillation around such a vev is typically both abundant 
■ and long lived, which modifies the simple picture in a significant and sometimes disasterous way. Extending an old 
On ' terminology [|l|, we shall call a scalar field with a large vev and a flat potential a 'flaton field', or simply a 'flaton'.|^ 
Although flaton fields are by no means inevitable, they are natural in the context of modern particle theory and in 
O ^ our opinion their possible cosmological consequences should be taken very seriously. Some aspects of the cosmology 
' . of flaton fields are already well known [l|J^-[To[| , and in a recent note ||l^ we drew attention to a new feature which we 
termed 'thermal inflation'. The present paper, along with two more in preparation ^2 13 1, aims to give a systematic 
(— I ' account of the subject. 

• • ■ Let us begin by being more precise about what is meant by a 'large' vev, and a potential which is 'almost flat'. 
. !^ ' These terms are defined with respect to the energy scale 10^ to 10"^ GeV, which is the scale of supersymmetry breaking 
! as defined by the masses of the supersymmetric partners of known particles ||l4| . The vev is defined as the position 
of the minimum of the potential, and a 'large' vev M is one satisfying AI ^ l(r GeV. An 'almost flat' potential V is 
5^ , one whose curvature is of order 10^ to 10^ GeV (except near any points of inflexion) out to field values much 

bigger than 10"^ GeV, and if the field has a large vev this is supposed to be true out to at least the vev. For an almost 
flat potential the particle mass m is therefore of order 10^ to 10^ GeV. From now on we drop the qualifier 'almost', 
referring simply to a flat potential. 

The most widely discussed flaton candidates are the moduli occurring in superstring theory. The potential of 
a modulus is indeed flat, and if its vev is nonzero it is typically of order the Planck scale Mpi = (SttG)^^/^ = 
2.4 X 10^* GeV. A modulus with such a vevQ is known to be fatal to the standard cosmology since the corresponding 
particles are very abundant and do not decay before nucleosynthesis pjlQ-Sn]. As we shall see, the failure to decay 



o 



^Note the etymology. The term 'flaton' refers to the flat potential, not to inflation. Conversely, the familiar word 'inflaton' 
refers to the field which is slowly rolling during inflation. We shall also use the term 'flaton' to denote the particle species 
corresponding to a flaton fleld. 

■^A field with these properties occurred in the first example of a nonrenormalizable supersymmetry-breaking hidden sector, 
which contained a single complex field. It was called the Polonyi field, and the associated problem [g| was called the Polonyi 
problem. Most of what we say concerning the moduli applies to any species with these properties. 
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before nucleosynthesis is likely to persist for any flaton with a vev exceeding 10"'^ GeV, making all such flatons fatal 
to the standard cosmology [Q. 

Moduli are by no means the only flaton candidates. On the contrary, any field (in the observable sector) with a 
vev much bigger than 10"^ GeV is likely to have a flat potential, and so to be a flaton. The reason, as we discuss in 
detail below, is that it is natural to construct all available mass scales from just the two basic scales m and M-p\. 
Apart from the moduli, the most familiar examples of fields with nonzero vevs are those which are charged under a 
continuous symmetry, the vev then indicating a spontaneous breakdown of the symmetry. If the symmetry is local 
then the field is by definition a higgs field, and presumably the examples of this type occurring in nature (apart from 
the higgs fields breaking electroweak symmetry) are the higgs fields breaking the GUT symmetry, whose vevs are of 
order 10^^ GeV. Alternatively the symmetry could be global, a likely candidate for this case being the Peccei-Quinn 
field with a vev perhaps of order 10^^ GeV. On the other hand, it makes perfect sense for a field to have a nonzero 
vev even if it is not charged under any continuous symmetry. For example, a right-handed neutrino mass might be 
generated by a vev, without lepton number being a good symmetry . 

As mentioned already, moduli as well as any other flatons with a vev bigger than 10^'' GeV are fatal to the standard 
cosmology. How are we to solve this 'moduli problem' if it exists? 

The usual recipe for getting rid of unwanted relics in cosmology is to invoke an e arly ep och of inflation, lasting 
at least 50 to 60 Hubble times or so. Such an era is also desirable for other reasons [ P2|j24| , one of which is that it 
can generate an adiabatic density perturbation of the right magnitude to explain the cosmic microwave background 
anisotropy and large scale structure. To do this the potential at the end of inflation must satisfy y^/** < 10^® GeV 
p4[ , and the lowest value of V^^'^ that has been proposed in a plausible model is V^/^ - lO^^ GeV 

Inflation at such a high scale does not solve the moduli problem, because although it sufficiently dilutes moduli 
present before inflation they are regenerated with an unacceptable abundance afterwards. We show in |[Tl|], and in 
much more detail below, that to avoid excessive regeneration one requires 



where Tr is the reheat temperature. An era of inflation at such a low energy scale seems impossible to realize in the 
context of sensible particle physics, if it is required also to produce the cosmological density perturbation. Randall 
and Thomas [l^ therefore suggested that the density perturbation is produced by an era of inflation at the usual high 
energy scale, while a second era of inflation at a low energy scale solves the moduli problem. However, even without 
the constraint of producing the density perturbation it is difficult to construct a model of inflation giving a sufficiently 
low energy scale, within the usual paradigm where there is an inflaton field rolling slowly down the potential. The 
reason stems from the fact that a necessary condition for slow roll is that the inflaton mass (or more precisely the 
curvature \V"\^/'^ evaluated while the field is rolling) be much less than the Hubble parameter H ~ l^^/^/Mpi. The 
bound displayed in Eq. (|l]) corresponds to a very low mass ^ 10 MeV. 

The central purpose of this paper is to explore the fact that a flaton fleld can lead to a completely different type of 
inflation, called thermal inflation JTH , which can solve the moduli problem provided that the vev M is within one or two 
orders of magnitude of 10^^ GeV. During thermal inflation the flaton field is held at the origin by finite temperature 
effects so that no field is rolling. The potential during thermal inflation is the value Vq of the flaton potential at the 
origin, which is of order m'^M'^. With M ~ 10^^ GeV this gives Vo^/" ~ 10^ GeV which can satisfy Eq. (0). Thermal 

inflation starts when the thermal energy density falls below Vb which corresponds to a temperature roughly Vq^^, 
and it ends when the finite temperature becomes ineffective at a temperature of order m, so the number of e-folds is 
i ln(M/m) ~ 10. It turns out that this can sufficiently dilute the moduli existing before thermal inflation (especially 
if reheating after thermal inflation is delayed) and it will not interfere with the density perturbation generated during 
ordinary inflation. There is also the intriguing possibility that two or more bouts of thermal inflation can occur in 
quick succession, allowing an even more efficient solution of the moduli problem. 

The present paper and its two successors are complementary to recent papers by Dine, Randall and Thomas 
1^,^. The latter focus on fields with a flat potential but zero vev. These fields too are liable to be oscillating in 
the early Universe and if they carry nonzero lepton or baryon number they can lead to baryogenesis (the Affleck-Dine 
mechanism). However baryogenesis in this way works only if there is no thermal inflation, and that in turn is a viable 
possibility only if there is no moduli problem. The two sets of papers therefore represent mutually exclusive scenarios 
for the early Universe, and only time will tell which if either is correct. 

The rest of this paper is divided into two main sections plus a concluding one. In Section 2 we study the effective 
potential expected for flatons, both in the early Universe and in the present era when it reduces to the ordinary low 
energy effective potential. Special attention is paid to the case of moduli, which is different from that of other flatons 
because the moduli potential vanishes if supersymmetry is unbroken. The flaton decay rate is also estimated. The 
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reheat process for homogeneous flaton oscillations is considered, taking account of possible parametric resonance. In 
Section 3 a systematic account is given of the history of the Universe, assuming that thermal inflation occurs and that 
there is a moduli problem. The concluding section summarizes the results, and points to future directions of research. 

II. FLAT POTENTIALS AND FLATONS 

In a generic supersymmetric gauge theory there will be a large number of directions in the space of the complex 
scalar fields^ in which the potential V is exactly flat, before supersymmetry breaking and non-renormalizable terms 
are taken into account. (This is true, for example, in the Minimal Supersymmetric Standard Model.) After these 
effects are taken into account the potential is still almost flat, in the sense that the energy scale \V"\^^'^ specifying 
the curvature of the potential is only of order 10^ to 10"^ GeV, out to field values many orders of magnitude bigger 
than this scale. In this paper we are interested in flaton fields, which by definition correspond to flat directions with a 
nonzero vev. The central theme of this paper is that flaton fields arc cosmologically significant, because they typically 
lead to thermal inflation, and because they in any case oscillate homogeneously until a relatively late epoch. 

A field with a nonzero vev is by definition either a higgs field or a gauge singlet. We will focus on the latter case 
in this paper, since a straightforward interpretation of the data indicate that the vev of the GUT higgs field is of 
order 10^^ GeV which is too high to give viable thermal inflation. Note, though, that in some GUT models there are 
additional higgs fields with much smaller vevs [2^ ]. 

The cosmology of a given flaton field is largely determined by the form of its effective potential. One needs to know 
both the low energy effective potential which is relevant at the present era, and the effective potential in the early 
Universe. Also, since the case of moduli is somewhat different from that of flatons in general we treat the moduli in 
a separate subsection after the general discussion. 

A. The low energy effective potential 

Consider a complex flaton field (j>. In the limit where the potential is absolutely flat there is a global U{1) symmetry 
under the transformation (f> — > e*"(/), with an arbitrary choice for the origin of (f). In the full theory this symmetry 
may survive for one choice of the origin, at least to a good approximation, or it may be so badly broken as to be 
unrecognizable. 

Global 17(1) symmetry 

We begin by considering the case where the symmetry survives. Extensions of the Standard Model can indeed 
contain spontaneously broken global Ujl) symmetries, a well known example being the Peccei-Quinn symmetry 
associated with the axion |2|,|2|,||,|^,|ll. We initially suppose that the U{1) symmetry is exact. The potential 
then depends on only through |0|, and assuming an effective theory that is valid right up to the Planck scale, the 
potential in the flat direction is typically of the form 

oo 

V = Vo- ml\q^\' + ^ A„Mpj2"|,/,|2»+4 (2) 

71=1 

The |0p term comes from soft supersymmetry breaking, which means that mg ~ 10^ to 10"^ GeV, and the higher 
order terms are non-renormalizable terms. The dimensionless couplings A„ are at most of order 1, if the theory is 
indeed valid up to the Planck scale. 

The crucial feature of this potential, which distinguishes it from the potential of a generic field and makes it flat, is 
the absence of a term Al^l"^ with A ~ 1. Such a term can be forbidden by discrete or continuous gauge symmetries, in 
combination with supersymmetry. Supersymmetry breaking then generates a X\(f>\'^ term with a suppressed coupling 



^Each scalar field is complex in supersymmetric theories because supersymmetry relates it to the two degrees of freedom 
associated with a left- or right-handed spin-half field. In this paper we are assuming that the fields are canonically normalized 
in the regime of interest. (One cannot in general canonically normalize the fields exactly over an extended region of field space.) 
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A ~ (mo/Mpi)^. Such a term is negligible for flaton fields which are not moduli and we have lost nothing by omitting 
it from Eq. (||). (The case of moduli will be discussed in a moment, and in more detail in Section 2.5.) 

As the notation suggests, we have in mind the case where the mass-squared at the origin, — TOq, is negative. This 
means that the vev of \(f>\ does not vanish but rather has a value M ^ mg. To estimate M, suppose first that all of the 
A's are of the same order. Then as one increases |(/)|, the term comes in first, leading to M = {3Xi)~^^^mQ^^ Mp^^ ^ 
X 10"'^'' to 10^"'^ GeV. Now suppose instead that this term is negligible, so that the term comes in first. Then 
M = (4A2)-i/6TOy^M^/' ^ A2-1/' X 1013 GcV. If more terms are absent the vev will be raised further so the predicted 
range is M > lO^" GeV. In the entire regime |0| ^ M the curvature |T/^"|i/^ of the potential is only of order mo, 
which is of order 10^ to 10^ GeV. In particular the mass m of the flaton particle is of this order, and from now on 
we shall generally use it instead of mo when writing down order of magnitude estimates. The requirement V{M) = 
gives Vq ~ m?M'^ , corresponding to 



\ 106 GeV / V IQi" GeV 
If the nth term dominates in Eq. (Q), thenQ 



m 



= 2(n+l)m^ (4) 
M2"+2Mpi2" = [2(71 + l)(n + 2)A„]- V2 (5) 
Vo^[2{7i + '2)]-'^m^M^ (6) 

Rather than the non-renormalizable terms being suppressed by the Planck scale, they might be generated by 
integrating out particles with GUT scale masses and so instead be suppressed by Mqut — 2 x 10^^ GeV. This would 
correspond to taking A„ < (-^^pi/-M^gut)^" and would give the somewhat looser lower bound 

M > 10^ GeV (7) 

We noted a moment ago that in Eq. (||) the term has a coupling A ~ (m/Afpi)^ which is many orders of 
magnitude less than 1. It may happen that the same is true of one or more further terms. But for a flaton which is 
not a modulus one expects to find, at not too high order, a term whose coupling A„ is not many orders of magnitude 
less than 1. As a result, one expects the vev of flaton which is not a modulus to be several orders of magnitude 
below Mpi. By contrast one expects for a modulus that all couplings are strongly suppressed, because the potential 
of a modulus vanishes exactly when supersymmetry is unbroken. A natural order of magnitude for the couplings of a 
modulus is (m/Mpi)^ making the vev of order Mpi, though there are other possibilities. We shall discuss moduli in 
more detail in Section 2.5. 

The flat potential Eq. (||) is not at all what cosmologists generally assume when they consider spontaneous 
symmetry-breaking in the early Universe. Rather they assume, as for instance in the textbooks ||2^,^,^ and the 
reviews that the potential is like the Standard Model higgs' potential^ 

V = A(|0|2 - m2)2 (8) 

with A ~ 1. For the Standard Model higgs, whose vev is of order 10^ GeV, this potential is indeed natural from the 
viewpoint of supergravity; it simply corresponds to a non-flat direction, in which there is a AJ</>|'* term. But when M 
is much bigger than m it becomes far less natural, and in our view Eq. (j^) rather than Eq. (|) should be regarded as 
the 'default' case. 

So far we have taken the U{1) symmetry to be exact, so that the goldstone boson corresponding to the angular 
direction is massless. If the symmetry is broken the goldstone boson will acquire a mass. This mass is by definition 
much less than that of the flaton if the symmetry is only slightly broken. On the other hand, as we now discuss 
the symmetry may be strongly broken which means that the would-be goldstone boson becomes just another flaton 
particle. 



*The mass-squared of the flaton particle is ^V"{M) because the canonically normalized complex field j^j is related to the 
canonically normalized real flaton particle field 50 by 101 = M -|- 5(j)/y/2. 

^In the case of the Standard Model is a doublet and the symmetry is SU{2) but this is an irrelevant complication for our 
purpose. 
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No U(l) symmetry 



As a simple example, consider the superpotential W — (A/4Mpi) 
corresponding potential is of the form 



with A ~ 1. After supersymmetry breaking the 



V{4>)^Vo-ml\(t>\'^ + \AW + . 



dW 



dW 



Vo-ml\<j>\' 



PI 



mi 

Ml, 



(9) 



(10) 



with mo and the magnitudes of A, B and C all of order 10^ to 10'^ GeV. 

In this example U{1) has been broken down to (which leaves 0"* invariant), and there are four vacua each with 

the same vev \(f>\ = M |A|~^/^my^A/p|'^. In a given vacuum there are now two particles with mass 10^ to 10^ GeV; 
one of them is the one corresponding to the radial oscillation that we considered before, and the other is the would-be 
goldstone boson corresponding to the angular oscillation. We shall generally refer to them both as flatons. Note that 
in the regime |0| <C M the U{\) symmetry is approximately restored, since the term — mpl^p dominates. 

The Zi symmetry surviving in this example has ensured that there are no linear terms in the expansion of 4> about 
the origin, and this feature will become crucial when we consider the effective potential in the early Universe. Of 
course any Z^ symmetry will do for this purpose, and it does not need to be exact. 

In our discussion — mg has been taken to be negative. If it is positive the potential has a minimum at the origin. 
If this is also the position of the vev (ie., if it is the absolute minimum) then the field is not a flaton and does not 
concern us. It can however happen, as for instance in the model of |Q| , that the origin corresponds to a false vacuum, 
with higher order terms generating a large vev so that we are dealing with a flaton. Thermal inflation with such a 
flaton is viable only if tunneling to the true vacuum is rapid, which is typically not the case. 

For simplicity we shall from now on make frequent use of the notation appropriate to the case where there is a 
i7(l), writing the potential as a function only of \4>\ and using m to denote the mass of the flaton particle. 



B. The flaton decay rate 



There is a general expectation that a flaton particle corresponding to oscillations around a vev M will couple only 
weakly to particles with mass much less than M . In particular, one expects ||l|,§,||j8| |l^] that the flaton decay rate T 
is at most of order /M"^ . 

Consider first the decay into a pair of identical spin zero particles which correspond to a real field V'l with the 
renormalizable effective interaction A|0p?/'^. Setting |0| equal to its vev this interaction gives a contribution 2AM^ to 
the mass-squared to^ . Barring a precise cancellation^] it follows that 

X<-{^)" (11) 

where the right hand side is at most ^(m/M)'^ or the decay would be forbidden by energy conservation. Substituting 
= M + 6(t)/\/2, one finds that the flaton decay rate corresponding to this interaction is 

A2 /M^' 



^=8^[m) (12) 

Maximizing this expression subject to the constraint Eq. ( p] ) gives F ^ 10~^TO"^/Af^. 

This effective interaction with a coupling of order A ^ (m/M)^ is quite natural. For instance an interaction 
might give some field X a mass of order M, and then an interaction ip'^X'^ would generate it through the diagram 
with a single X loop. 



®It has been pointed out to us by G. Dvali that such a cancellation does occur in an SU(5) GUT where the doublet-triplet 
splitting problem is 'solved' by a fine tuned cancellation. In such a case the decay rate has the unsuppressed value F ~ m. 
When the problem is solved in a more acceptable way this need not be so, but we will not pursue the point here because our 
main focus is not on the GUT. 
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For an effective interaction involving more powers of the fields and/or derivatives the arguments are generally less 
precise, but one expects suppression because such terms are non-renormalizable and therefore involve inverse powers 
of some scale M which is presumably at least of order M . Consider for instance a term involving one power of and 
two of f/), with 2n derivatives. Its coefficient is expected to be at most of order X'M~^" with A' ^ M, and since the 
energy of all particles is of order m (in the (p rest frame) this gives the decay rate Eq. ( p^ with A ^ (m/M)^". For 
n > 1 this is much smaller than the upper limit Eq. (|l l[) , but for ?i = 1 it is bigger by a factor (m/m^)^ leading to 
r ^ (87r)^^r7i'^/M^. On the basis of this discussion, we shall assume that 

r = IQ-^-fm^/M'^ (13) 

with 7 1. 

The decay into goldstone bosons 

A definite example of a derivative coupling is provided by the decay of the 'radial' flaton into the 'angular' fiaton, 
or goldstone boson. Near the vev, the canonically normalized radial field s and angular field a are defined by 

= (^ + M) exp(m/V2M) (14) 
v2 

Expanding the canonical kinetic term £kin = df^<j}*d^(l) to first order in s, one finds the canonical kinetic terms for s 
and a plus an interaction term 

Ant = -^^d.ad^^a (15) 

The coefficient is of the advertised form A'/M^, with A' — M/^/2. 

The goldstone bosons produced by this coupling can be cosmologically dangerous, because their interaction can be 
too weak to thermalize them. This will be discussed in connection with the axion in (see also ]3l|]). 

The flaton freeze-out temperature 

Though we have focussed on the decay rate, similar considerations apply to collision rates. The rates for collisions 
involving a fiaton and other light particles are suppressed at energies well below Af , and therefore the freeze-out 
temperature below which fiaton particles cease to be in thermal equilibrium is very roughly of order M . Note that 
this applies only in the true vacuum, where the flaton field is oscillating about the vev. 

C. The effective potential in the early Universe 

In the early Universe, the interaction of a given field with other fields will alter the effective potential of that field, 
and in particular the effective flaton potential V{(j)) will be altered. 

We should flrst clarify what is meant by the 'effective potential T^(0)'. There is in reality a single effective potential 
V{4)^ ip, . . .), which is a function of all the scalar fields. It is natural to define the effective potential of any individual 
field as the full potential with all other fields held at their vevs, and this is the definition that we had in mind for the 
low energy effective potential V {(/)). However in the early Universe all sufficiently light scalar fields are significantly 
displaced from their vevs, either homogeneously in the manner we have been discussing for fiatons, or inhomogeneously 
as for instance if the field is in thermal equilibrium. Instead of evaluating the full effective potential V{(f>, ^, . . .) with 
the other fields at their vevs one should set them equal to their current time-averaged values, so that for instance 
a term is replaced by {ilP')<jP . In addition, the actual form of the full effective potential is affected by the 

presence of particles with nonzero spin and also by kinetic terms, so that the effective potential Vaaiyii'i 4'^ ■ ■ ■) in the 
early Universe is different from the low energy effective potential t4ow(<^, "0, ■ • •) which applies at present. For both of 
these reasons, the effective potential Vcai\y{(f>) in the early Universe is different from the low energy effective potential 
Mow(0) which applies at present. 

Although the form of the effective potential V{(l)) changes with the history of the Universe, its gradient will always 
vanish at the origin provided that it is invariant under at least a Z„ symmetry. This tends to be at least approximately 
true in simple models, and we shall take it for granted in what follows. Let us pause briefly though to see why such a 
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symmetry is common. If the full potential V{(j), tp, ■ ■ ■) is expanded as a power series in all of the fields each individual 
term will be invariant under one or more Zn symmetries unless it consists of just the first power of one field. For 
instance the term (jx^ijj'^ is invariant under a Z2 acting on (j), and another acting on ip. As we discussed in Section 
2.1, only a few leading terms will be important in practice, so it is reasonable that one or more Zn symmetries 
will be approximately present in the full potential. Then the question of whether or not the potential V{4>) of an 
individual field possesses an approximate Z„ symmetry depends on the form of the full potential, but again this is 
not unreasonable. 

Taking it for granted that the gradient of vanishes at the origin, let us ask what is the effective mass-squared 
y"(0) in the early Universe. (We continue to assume for simplicity that there is a U{1) symmetry, so that V is a 
function only of \(j)\.) 

First consider the era of ordinary inflation. It has been known for some time [pp,p|,^ that by looking at the form 
of the full potential predicted by TV = 1 supergravity one can identify contributions of order ±H^ to the mass-squared 
of every field. For the inflaton field(s) these contributions have to cancel because otherwise inflation will not occur, 
but for a generic field one does not expect a cancellation. Assuming that flatons are not inflatons, the conclusion is 
that their mass-squared during inflation is (at least) of order ±H^. 

After inflation it is not so clear what the mass-squared will be. In the extreme case where the interaction is of 
only gravitational strength one expects a contribution of the same order, We noted earlier that in the 

true vacuum, the interaction of flaton particles with other light particles is suppressed, so one at first sight expects 
something like this estimate to hold for a fiaton field. However, that suppression occurs because the vev of the flaton 
field is large (the flaton particles correspond to small oscillations around the vev) . Near the origin the flaton field can 
have unsuppressed interactions with light fields. 

To see why, take as an example the interaction iA|0pV^ that we considered earlier. When (j> is at its vev this gives 
a contribution AM^ to rrA. Barring cancellations, A must therefore be small if is small. But suppose that in 
contrast is of order M and is generated by this interaction. Then there is a coupling A ~ 1, and for flaton field 
values near the origin the -0 field becomes light. The result is that near the origin the flaton field has an unsuppressed 
interaction with the light field ^p. 

If is a higgs field, charged by definition under a gauge symmetry, a coupling of this kind to at least the gauge 
bosons and gauginos is inevitable. In the case where is neutral under all gauge symmetries, which is our focus here, 
such a coupling is not inevitable but it is still quite natural; for instance, in models of the kind discussed in |pl^,0,n2| 
a flaton field couples in this way to the right handed neutrino and sneutrino. 

Assuming that the flaton field near the origin indeed has unsuppressed interactions with one or more particle species 
having effective mass of order |0|, it will be in thermal equilibrium in the regime \(f)\ < T. (The upper limit comes 
from the fact that at a given temperature particles with mass bigger than T become too rare too maintain thermal 
equilibrium.) The finite temperature correction to the effective potential gives the flaton an effective mass-squared 
p],p5[ of order (T^ — TOq), which gives the effective potential a local minimum at the origin for T bigger than some 
critical temperature Tc mo ^ m. (As usual, —rriQ denotes the effective zero-temperature mass-squared at the 
origin, and m denotes the flaton particle mass which is the parameter we normally focus on. Recall that both toq and 
771 are of order 10^ to 10^ GeV.) 

In addition to the local minimum at the origin, the effective potential retains its true minimum at (f) — M except 
at very high temperatures T ^ Af , but there is no significant tunneling between the two [p],p5|. 

To summarize this discussion, if the fiaton field has gravitational strength interactions its mass-squared is expected 
to be of order ztH^. If, on the other hand, it has unsuppressed interactions then it will be in thermal equilibrium in 
the regime |0| < T and in this regime there will be a necessarily positive mass-squared of order ~ {Mp\/ H)H'^ . 
These are the most important possibilities for the effective mass-squared but others exist, especially during inflation 
where one might have a coupling to the inflaton field, say of the form fp'^'pfnf particular, hybrid inflation [p5| makes 
essential use of such a coupling). As in this example, the positivity of the potential tends to require that such a 
coupling again gives a positive mass-squared. 

D. The cosmology of fields with flat potentials 

In the light of what we have done so far there are the following four possibilities for the cosmology of a field with 
a flat potential. 

(i) The field sits at the origin. If the minimum of the potential is at the origin throughout the history of the 
Universe then the field will sit there apart from thermal and quantum fluctuations. In that case it does not undergo 
homogeneous oscillations in the early Universe, and we are not concerned with it here. It will in general have 
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unsuppressed interactions (at least if it is not a modulus) and the corresponding particle species will be produced 
through particle collisions and decays involving these interactions. 

(a) The field oscillates about the origin. Now suppose that although the minimum of the low energy effective 
potential is at the origin, the minimum in the early Universe is displaced because there is a negative mass-squared 
of order — i?^. In that case the field will start to oscillate about the origin at the epoch H ^ m. The oscillation is 
generally short lived, because the particles corresponding to it generally have unsuppressed couplings (except perhaps 
for moduli). If there is no thermal inflation the oscillation can however lead to viable baryogenesis through the 
Affleck-Dine mechanism [p7| . 

(Hi) Thermal inflation occurs. In the two remaining cases the vev is nonzero, so that we are dealing by definition 
with a fiaton field. Thermal inflation, which is the focus of the present paper, occurs if the fiaton field is held at zero 
in the early Universe by the finite temperature. It ends when the temperature falls to some critical value Tc ~ m 
(provided that the zero-temperature effective potential has no barrier separating the origin from the vev), after which 
the fiaton field starts to oscillate about the vev. The oscillation around the vev might persist for a long time because 
the coupling of fiaton particles to other light particles is suppressed (Section 2.2, and Section 2.6 below). 

(iv) Flatons not leading to thermal inflation. In the fourth case the flaton field fails to be held at the origin by 
the finite temperature of the early Universe. This will occur if the fiaton has an effective mass-squared ~ —H^ which 
prevents it from ever being near the origin. It will also occur whatever the sign of the mass-squared, if the interaction 
of the flaton field is suppressed even near the origin. When H falls to a value of order the flaton mass to, the held 
starts to oscillate about the vev, with an initial amplitude of order M . (The initial amplitude is equal to M if the 
initial field value is at the origin. If the field is displaced from the origin by a mass-squared of order —H^, its value 
is typically of order M when the field starts to oscillate.) As in the previous case the oscillation might last for a long 
time. 



What we have done so far, including the summary of the last subsection, applies in essence to all flaton fields 
including any which are moduli. On the other hand, moduli do have some properties which distinguish them from 
other scalar fields ('matter fields') and as a result the general discussion acquires a somewhat different flavour when 
applied to them. 

The low energy effective potential of a modulus vanishes exactly if supersymmetry is unbroken. After supersymmetry 
breaking its potential is generally thought to be flat, so that its curvature is everywhere of order 10^ to 10^ GeV 

(except near points of inflexion) . If a modulus has a nonzero vev, then as we discuss in a moment its vev is generally 
expected to be of order Mp\. To a large extent its properties can then be obtained simply by setting M = Mp\ 
in formulas that apply to flatons in general, but there are some special features. These arise because one is forced 
to consider field variations of order Mpi, in contrast with matter fields where one need only consider much smaller 
variations (typically of order the vev M <C Mpi for a fiaton field which is not a modulus). 

In order to talk about a nonzero vev for any field there has to be a well defined origin, which will be defined as a 
point which is invariant ('fixed') under the group of symmetries under which the field transforms. For matter fields 
this defines a unique origin, such that the symmetry group consists of linear operators in field space. For moduli the 
symmetries are more complicated, and there are in general an infinite number of fixed points with a separation of 
order Mp\ (though only a finite number are physically distinct because the symmetry is a discrete gauge symmetry) . 
If the vev of a modulus is at a fixed point it is natural to say that it vanishes, and otherwise it is natural to define the 
vev as the distance to the nearest fixed point. These are the conventions that we have had in mind, without explicitly 
stating them. The statement that the vev of some modulus is of order Mpi just means that it is not close to any 
particular fixed point. As with other fields, a modulus can have unsuppressed interactions with other light fields only 
if it is close to a fixed point. 

Each of the four possibilities for the cosmology of a flaton field listed in the last subsection exists for a modulus. If 
possibility (i) holds for all moduli then there is no moduli problem. Assuming that this is not the case, let us look at 
the expected form of the effective potential of a modulus $. For simplicity we will pretend that $ is real, and take 
it to be canonically normalized. Before supersymmetry breaking is taken into account the potential V{^) vanishes. 
With the breaking taken into account the potential in the true vacuum (the low energy potential) is generally thought 
to be of the form 



Here the supersymmetry breaking scale Ms is related to the scale m ^ 10^ to 10'^ GeV by Ms ~ (mMpi)^/^ ^ lO^*' to 
10^^ GeV, and f{x) is a function whose value and low order derivatives are typically of order 1 in the regime \x\ < I. 



E. The moduli potential 




(16) 
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We have expanded the potential about its vev $i. Note that the potential vanishes in the limit Ms of unbroken 
supersymmetry, in accordance with the fact that we are dealing with a modulus. 

In the early Universe there will be additional supersymmetry breaking because of the nonzero energy density p, 
leading to an additional contribution to the potential of the form P,^, ^ , po| 

Kosm-p.9(^^) =^i?'(<i>-<i>2)' + ... (17) 

The function g{x) has value and low order derivatives of order 1 (making a ^ 1). The minimum of this potential is 
located at a different value $2, which is displaced from the true vev $1 by a distance ^0 = ^2 — ^1 ^ Mp\. 

We have in mind the case where both $1 and $2 are nonzero (case (iv) of the last subsection) and of order Mp\ . 
If $2 = but $1 7^ (case (iii)) there is also a moduli problem, but it might be rendered insoluble by domain walls 
(though in analyzing this possibility within a given model one will have to remember that the discrete symmetries 
under which the moduli transform are gauge symmetries). If $1 =0 but $2 7^ (case (ii)) there is no moduli problem 
if the relevant moduli have unsuppressed couplings near the origin. 

Although Eq. (^6|) is the simplest possibility for the potential of a modulus there are others, which could lead to 
a vev below the Planck scale. For example, if supersymmetry breaking is due to hidden sector gaugino condensation 
then the moduli potential might include terms of the form ^"|(/)|'"+^/Mp[^'" where ^ is related to vevs arising from 
gaugino condensation. These terms still vanish when supersymmetry is unbroken, as is required for a modulus, but 
they might generate a vev below the Planck scale. For example the GUT Higgs could be a Wilson line modulus, with 
a vev of order 10^^ GeV generated in this way In considering the moduli problem we assume in this paper that 
at least some moduli have a vev of order Mpi. 



F. The flaton reheat temperature 

Let us quantify the statement that the flaton field oscillations in the early Universe last for a long time. 

The oscillation of a flaton field with vev M has initial amplitude (po ^ M. The corresponding energy density is 
Pij, ^ ^TO^^Q, and the number density of the flaton particles is ~ ^m(j)Q. These particles have no random motion 
because the field is homogeneous, so they constitute matter as opposed to radiation. If the flaton is associated with 
thermal inflation, the oscillation commences after thermal inflation and immediately dominates the energy density. 
If not, the oscillation commences at the earlier epoch H ^ m, and may or may not come to dominate the energy 
density. 

If the oscillation amplitude decreased like a^^^^, where a is the scale factor of the Universe, then the energy per 
comoving volume of the flaton field would be conserved. In fact, the energy drains away through the interactions of 
the flaton field so that the oscillation amplitude decreases faster. 

If the oscillation amplitude is sufficiently small and the interactions are sufficiently weak, each fiaton particle decays 
independently so that the rate at which the energy drains away is simply the particle decay rate P. It has practically 
all disappeared soon after the time 

T. 1 o i.r. 9f M \V300GeV\^ 

"^^'^ [wh^) [-^) '''' 

where we have used Eq. (p3[). Setting this time equal to H^^ and assuming that the decay products thermalize 
promptly we arrive at an estimate of the 'reheat temperature', 

lO^GeV^ ' - ^ ' 



To . ,-PtM|, . 3,t ( ' GeV (19) 



where g^, ^ 10^ is the effective number of species at T = TdJ^ 



''The following results will be used without comment in the text. The entropy density of radiation at temperature T is 
s = (4/3)p/r = (27r^/45)g,T^ = 1.0X3^''^^''*, where g,(T) is the effective number of particle species in thermal equilibrium, 
and p = (7r^/30)(ji»r'* is the energy density. As the Universe expands the scale factor a increases. The energy density in 
relativistic particles (radiation) is proportional to a""* and that in non-relativistic particles is proportional to a~^. In thermal 
equilibrium the entropy a'^s in a comoving volume is constant and so is gV^oT. According to the Standard Model, gV^ is in the 
range 1 to 2 for T < 100 MeV, and then rises sharply to become ~ 3, finally rising to ~ 4 when T > 10"^ GeV in supersymmetric 
extensions of the Standard Model. We use the appropriate value in our estimates. 
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As has been discussed recently in connection with ordinary inflation, the assumption that each flaton particle decays 
independently need not be correct ||37|-^ (see also Instead, parametric resonance effects can drain away much 

of the oscillation energy as soon as the oscillation starts, leaving behind only some fraction to decay at the single 
particle decay rate. The energy drained away goes initially into the creation of marginally relativistic scalar particles. 
All species are produced which have sufficient coupling to the flaton, including the flaton itself. (We are not aware 
of any discussion of the possibility of the production of bosons with spin 1 or higher through parametric resonance 
and it may be that this also occurs. Fermions are not produced in signiflcant number because of Pauli blocking.) If 
nothing happens to the produced scalar particles they will become non-relativistic after a few Hubble times, and are 
expected to decay at their one-particle decay rate.|^ If, on the other hand, they thermalize then they turn into highly 
relativistic radiation. 

At the present time it is not clear whether parametric resonance can really create particles which thermalize 
successfully. However, it is clear that the flaton component of the produced particles cannot thermalize because here 
one knows that the interaction is too weak. Furthermore, one expects that the energy density of the produced flatons 
will be a significant fraction of the total energy density |Q. Thus, even if the other produced particles thermalize 
promptly one expects that a significant fraction of non-thermalized energy will remain, and that a significant fraction 
of that energy will be in flaton particles. 

Any thermalized radiation produced by parametric resonance will redshift away, so independently of the details one 
expects that a few Hubble times after the end of thermal inflation the energy density is dominated by non-relativistic 
scalar particles, including the flatons and perhaps other species. Each species will decay at the single-particle decay 
rate, so we expect eventually to find only the longest-lived species, which dominates the energy density until it decays. 

For simplicity we shall assume in what follows that this species is the fiaton itself, and we shall also ignore the effect 
of any radiation produced by particle decay. Thus we are in effect assuming that soon after thermal infiation has 
ended, some fraction e of the energy is in non-relativistic flaton particles which decay according to the one-particle 
decay rate, with the remainder in thermalized radiation. This should describe the real situation at least approximately, 
provided that any non-flaton particles produced decay at least as rapidly as the flatons. The important special case 
e = 1 is considered, and the possibility that e may be very small is not discounted. However, as as discussed above, 
this latter case seems unlikely because one expects that parametric resonance will convert a significant fraction of the 
energy density into flaton particles which interact too weakly to thermalize. 

The upshot of this discussion is that despite the possible occurrence of parametric resonance, one expects that 
the eventual reheat temperature after thermal inflation is still the temperature Tq calculated from the single-particle 
decay rate, as given by Eq. (|l^). If Td is indeed the reheat temperature, the requirement that it be not too low 
places strong restrictions on M. In order not to upset nucleosynthesis one must have Td ^ 10 MeV, which requires 
M < lO^^'GeV (taking m < lO^GeV)]^ However, if R parity is respected as is usually supposed, there is a stable 
LSP which imposes a much stronger constraint. Indeed, to bring the LSP into thermal equilibrium so that it is not 
over-produced (and can naturally have the correct abundance to be the dark matter), one needs Td substantially in 
excess of the LSP decoupling temperature which is of order 1 GeV. Thus one needs M < 10^^ GeV. Finally, one might 
wish to generate baryon number through the electroweak transition which would require Td ^ 100 GeV corresponding 
to M ^ 10^" GeV. In view of the fact that these limits are perhaps rather conservative (since one expects 7 to be 
significantly less than 1, and does not anticipate m as high as 10^ GeV) this last requirement is hardly likely to be 
satisfied, but other baryogenesis mechanisms exist as discussed in [Oj. 



When the the particles have become non-relativistic one might think that parametric resonance will recommence, since 
the wavenumber of the corresponding scalar field is negligible compared with its frequency. However, the collection of non- 
relativistic particles corresponds to a superposition of almost-classical quantum states, not to any one such state, since the 
phases of the corresponding fields are uncorrelated, so it is not clear that the parametric resonance formalism applies. More 
importantly, the amplitude of the would-be classical oscillation will typically be too small for parametric resonance to occur. 
We are indebted to A. D. Linde for helpful correspondence about this issue. 

^In we estimated M < lO^'^ GeV. The extra factor 100 came from three different sources. First we used the very naive 
estimate F ~ rr? jM^ , corresponding to 7^''^ = 10. Second, we set — 1 where as the true value is more like 10^/^ Third, 
we rounded up our estimate of Td to the nearest power of ten which meant multiplying it by of order lO^''^. It so happened 
that each of these approximations went the same way to give the factor 100. 
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III. COSMOLOGY WITH THERMAL INFLATION 



We now give a systematic account of the history of the early Universe in the case where there is thermal inflation. 
We assume that there is a moduli problem because this provides the strongest motivation for thermal inflation, and 
assume that at least some of the moduli have a vev of order Mpj. We also assume that any radiation produced 
by parametric resonance promptly thermalizes. With these assumptions there are the following eras which we shall 
consider in turn. 

1. Ordinary inflation. 

2. Matter domination by the homogeneous oscillation of the inflaton (unless full reheating occurs promptly). 

3. Full reheating, which leads to radiation domination if it occurs before the moduli start to oscillate. 

4. Homogeneous oscillation of the moduli, starting at the epoch H ~ r7i$. If reheating has previously occurred 
there is now matter domination by the moduli. If it has not occurred the moduli and inflaton matter densities are 
roughly comparable, and remain so until full reheating (of the inflaton matter). We assume that full reheating 
takes place before the beginning of thermal inflation. 

5. Thermal inflation. 

6. Matter domination by the homogeneous oscillation of the flaton field which caused thermal inflation (unless 
reheating occurs promptly). 

7. Full reheating of the flaton matter, leading to radiation domination before nucleosynthesis after which the history 
of the Universe is the standard one. 



A. Before thermal inflation 



One expects the Universe to start with an era of ordinary inflation [£2|,g3| , whether or not there is a later epoch of 
thermal inflation. During this era, the energy density p is dominated by the potential V of the scalar fields, with all 
except the inflaton field (or fields) fixed. The infiaton field slowly rolls down the potential, because in its direction 
the flatness conditions |MpiV^'/T^| <^ 1 and \V"\ <^ are satisfied 22, 2^]. We noted earlier that in the context of 



supergravity the second of these conditions requires cancellations. Although these might be accidental it is attractive 
to suppose that they occur by virtue of some symmetry. One suitable symmetry (most easily implemented in the 
context of hybrid inflation |^^) was suggested in and another has been proposed in |Q. A third possibility 

is to invoke a global U{1) symmetry as in but this is problematical because the inflaton potential vanishes in the 
limit where the symmetry is exact so that the magnitude of V" is difficult to control]^ 

To avoid generating too much large scale cmb anisotropy the potential at the end of ordinary inflation must satisfy 
11 

V^''^ < 10^^ GeV (20) 

At some epoch after ordinary inflation 'reheating' occurs, which by definition means that practically all of the 
energy density thermalizes (except for the contribution of moduli). If reheating is prompt the reheat temperature is 
7r ~ {y/g*y^^- A naive estimate of the time taken for reheat would be that it is the decay time of a single infiaton 
particle, which typically leads to a much lower reheat temperature. However prompt conversion of a large fraction 
of the energy density into marginally relativistic particles is likely. In the commonly discussed case where inflation 
ends with the oscillation of a homogeneous inflaton field this is expected to occur through the parametric resonance 
effect that we considered already for the case of thermal infiation. It is also expected to occur in the case of hybrid 
inflation though a quantitative account of this case has not yet been given, and will be more complicated because 



^"An alternative idea is to suppose that the potential is exactly flat (or at least much flatter than that of the inflaton 
field) in the direction of at least one field, say a modulus, which couples to the infiaton. The inflaton potential then depends 
on the value of this field, which will vary from place to place in the Universe allowing the possibility that we live in a region 
where the infiaton potential happens to be sufficiently flat. But this just pushes back to another level the problem of finding 
cancellations which keep the potential flat in some direction. 
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spatial gradients are probably important from the beginning pGf . As we discussed earlier these marginally relativistic 
particles may then thermalize promptly leading to full or partial reheating .p] 

Moduli (more precisely, those moduli if any which are flatons with M ^ Afpi) are produced both before and after 
thermal inflation, and we shall call the moduli from these sources respectively big bang moduli and thermal inflation 
moduli. 

When H ^ m$ the modulus' potential is given by Eq. ([l7|), so that $ is shifted from its true vacuum value by 
$0 — ^2 — ^1 ~ Afpi. $2 will depend on the composition of the Universe and so $0 will change at any phase 
transitions, such as the end of inflation, but $ will rapidly settle down to its new minimum as it is critically damped. 
However, at the epoch H ^ m$ it starts to oscillate about the minimum of its low energy effective potential, and 
after H has fallen significantly below to$, the oscillations will no longer be critically damped and so are much more 
dangerous. 

During thermal inflation H <C m$, so the effective potential is dominated by Vtruc but ycosm still gives a small 
contribution, so the position of the minimum is shifted slightly from the true vacuum value. Oversimplifying a bit, 
we can estimate the shift by adding together Vtruo and Vcosm which gives 

y = im|(<i>-$i)2 + |i?2(<i>-<i>2)' + ... (21) 

= ^mld<f' + ^H^6<i>-<Pof + --- (22) 

= i (ml + aH') (S<P - , $0) ' + . . . (23) 

where (5$ = $ — $1 is the displacement of $ from its vev. In the last line a is of order 1, so the minimum of the 
modulus' potential is shifted during thermal inflation by an amount of order {H/m^)^Mpi pH |. 

To estimate roughly the abundance of big bang moduli, we can assume that the modulus field starts to oscillate 
about its vev when H ^ m$ with amplitude of order $0 ~ Mpi. The energy density ~ m|$Q/2 is of order the 
total energy density. If reheating has already occurred one can crudely set the radiation energy density equal to the 
total energy density which leads to the estimate 

s lOMpj m^/^ 

(In this expression s is the entropy density, and we are using the standard results summarized in the footnote after 
Eq. (^9|).) If reheating occurs later the moduli energy density is a fixed fraction of the total until reheating, and again 
setting the radiation density equal to the total density after reheating one finds 



2^1/2 

(25) 



10Mp(V<i> 



It is described in the Appendix how a more sophisticated calculation leads to the same results. 

We shall assume that full reheating occurs before the onset of thermal inflation (except for the contribution of 
moduli). The opposite case will be discussed in [p^ . 

These estimates for the moduli apply to any flaton not giving rise to thermal inflation (option (iv) of Section 2.4), 
if $0 is replaced by M. 



B. Thermal Inflation 



Thermal inflation will occur if one or more of the flaton fields is trapped at the origin in the early Universe. For 
the moment we suppose that only one is trapped. 



^'^Note, though, that an extremely low fraction of the energy density cannot thermalize because thermalization requires that 
the interaction rate per particle exceeds H. If the decay products are charged under some gauge symmetry, this requires 
ct9*T ^ H where a is the gauge coupling. Setting a ~ g*^^^ ~ 10~^, one finds that it is satisfied only if the fraction is bigger 
than V/{1Q^^ GeV)*. This constraint does not seem to have been noted before in the literature. 
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The trapping may initially be due to a non-thermal contribution to the mass-squared such as that of order H^. 
However if full reheating occurs before the beginning of thermal inflation then well within a Hubble time of the end 
of inflation enough entropy to trap the flaton at zero will have been released even by the single particle decay of the 
inflaton. 

If full reheating is indeed delayed to the epoch when thermal inflation begins, the temperature at that epoch is of 
order ^^'^Vq^^ ^ (jnAI)^^'^ corresponding to 

V 106 GeV y VlOi"GcV 

At the other extreme where reheating occurs before the moduli start to oscillate, the temperature at the beginning 
of thermal inflation is reduced by a factor (M/Mpi)^/6_ During thermal inflation T oc exp(— iJi) and it ends at 
T = Tc ^ 771, so there are at most of order ^ ln(M/m) ~ 10 e- folds of thermal inflation. This will not much affect the 
cosmological density perturbation generated about 50 e-folds before the beginning of ordinary inflation, though there 
might be a slight change in the spectral index. 




C. Entropy production after thermal inflation 

After thermal inflation ends, relic radiation from the first hot big bang plays no further role. The flaton field now 
starts to oscillate around its vev with initial amplitude M, corresponding to non-relativistic flatons (matter) which 
dominate the energy density. 

The decay of the flaton field generates entropy. If there is no parametric resonance the entropy per comoving volume 
increases linearly from the end of thermal inflation until the flaton decays, leading to an increase in the entropy by a 
factor 



(27r2/45)5,(Tc)T3 TSTdT^ 

Now suppose instead that there is parametric resonance which promptly thermalizes a substantial fraction of the 
energy density, leaving a fraction e in the flatons. This will increase the entropy by a factor 



(2^2/45)5* (Tc)T3 25T^ 
10 / M \^ f fSOOGeYY f Vq 



(28) 

Uoi"Gcvj [jh) [ m ) ['^^i^) ^^^^ 

The radiation energy density may initially dominate, but we assume that it falls below that of the residual flatons 
before the epoch To. The entropy release from the decay of these flatons is signiflcant only during the era Td ^ 



T < (TpTeq)^/^ 1 23 1, so it is a good approximation to regard this entropy release as suddenly occurring at the epoch 
T ^ Ty). It increases the entropy by a further factor 



(30) 



106e ( ] " (9^] ( 1 ^ ) (31) 

VlO^"GeVy \Tu J V300GeV/ ym-^APj ^ ' 



The total entropy increase is 



A ~ Apr Ad ^ (32) 



c 

)i5-5e ( ^'^ V /GeV\ / 300 GeV \ / Vn 



Eq. (^) is only supposed to apply if it gives a value Ad bigger than 1, which fails to be true in the small e 
regime e < Tu/Vq^^. This is the regime in which the flaton oscillation fails to dominate the energy density before it 
disappears at the epoch T = Td- In it Ad is practically equal to 1, and A = Apr has the e-independent value given 
by Eq. (|8l). 

We shall not consider the case where parametric resonance creates radiation which fails to thermalize, and hence 
quickly reverts to matter in the form of homogeneously oscillating scalar fields. 
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D. Solving the moduli problem with single thermal inflation 



In order not to upset nucleosynthesis, the moduli abundance n^/s must be less than 10 to 10 when nucle- 



osynthesis begins 49 1. Let us see what is required to satisfy this bound, first for the big bang moduli and then for 
the moduli produced after thermal inflation. 

We can assume that the flaton oscillation comes to dominate the energy density, because the assumption can be 
shown to be valid in the regime of parameter space satisfying the nucleosynthesis bound on the moduli abundance 



and to lead to an overestimate of the moduli abundance outside this regime. As a result we can use Eq. (32), and 



combining it with Eq. (25) one finds that the abundance of big bang moduli after thermal inflation is 

(34) 



10-16 



IQi^ GeV y /1\ / MpiHn y ( Tp \ 



m^J [mpJ VSOOGcV/ V Vq 



(35) 



In these formulas Hp is to be considered as being in the range m$(M/Mpi) ^ < m$. The lower limit comes from 
our assumption that full reheating after ordinary inflation occurs before the beginning of thermal inflation, and if Hp 
actually exceeds the upper limit the above formulas give the correct result when it is set equal to this limit. 

To analyze these constraints, assume first that lb > 1 GeV as is required if the LSP is stable, and recall that this 
implies M < 10^^ GeV, from Eq. (pj|). In Eq. (p5|), the round brackets in the second line are all of order unity, so we 
see that the big bang moduli may be sufficiently diluted for M as low as 10^ GeV, though this requires all parameters 
to be pushed to the limit and a more reasonable estimate of the lower limit might be 10^"'^ GeV. Now assume only 
that Td ^ 10 MeV, as required by nucleosynthesis, which implies M < 10^'' GeV. Then we see that unless e is very 
small it should be possible to solve the moduli problem, with no significant additional constraint on M. 

Now consider the moduli produced after thermal infiation. From Eq. (p3|), the minimum of the potential during 
thermal inflation is displaced from its true vacuum minimum by an amount (5$ ^ (Vo/"^|,Mpj)$o- The dynamics 
at the end of thermal infiation will be complicated but one would expect to generate a moduli number density 
n<i> ^ m$(5$^/2 ^ ^qVq /2m^Mp^ at the end of thermal infiation. Therefore the abundance of thermal infiation 
moduli is expected to be 

n$ _ 'PlVi/2mlM^, _ 'fgFoTD ^^^^ 



s 



10 



-15.5 



M \ V 1 \ / Td 



1012 GeV / \e \ GeV 



$0 V /300GcV\ / Vo 



Mp\ I \ m$ 



(37) 



Bearing in mind the relation between Td and M, we see that the abundance of thermal moduli does not impose a 
significant additional constraint. 

Moduli will be produced in the fiaton's decay with abundance 

s V s 

Since the flaton energy density is toti^ and we are assuming that it all thermalizes, rt^/s is of order Td /m and therefore 

s ^ \Tu J V"^V87^M2J V 300 GeV 



^0-'" oZlr. (39) 



which is probably sufficiently small. 

Finally we consider the possible thermal creation of gravitinos, moduli and modulinos after thermal inflation. 
Gravitinos, for which the most detailed calculations exist, appear to be created in a cosmologically safe abundance 
provided that the maximum temperature is less than | |50[ | 10^ GeV, and a similar result presumably holds for moduli 
and modulini since in all cases the interaction with other particles is of gravitational strength. This bound is satisfied 
after thermal infiation even in the extreme case where most of the energy density thermalizes immediately. 
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E. Double Thermal Inflation 



So far we assumed that only one flaton field gives thermal inflation, or in other words that only one flaton field has 
a thermal mass-squared which traps it at the origin in the early Universe. If two or more flaton fields are trapped the 
situation is in general much more complicated, but it simplifies considerably if the fields do not interact significantly. 
We treat this simple situation now, leaving the case of interacting fields to future publications jl^Jl^. Thus we 
consider two flaton fields and (j)2 , and assume that each of their potentials is of the form Eq. @ , 

y(0i,02) = Vi+V2~ ml\cj)i\^ - ml\^2\^ (40) 

The higher order terms stabilize the fields at ipi — Mi, and the constants Vi and V2 are the values of the separate 
potentials at the origin, with Vi ~ mlMf . The critical temperatures at which the fields roll away from zero are Tci, 
and we take Tci > Tc2- When the temperature drops below Tci, (pi will roll away from zero. 

If parametric resonance does not produce significant thermalization, the second field now also rolls away promptly 
and the situation is not substantially different from the case of thermal infiation. If on the other hand a significant 
fraction of the energy density is thermalized by parametric resonance, the temperature will be raised sufficiently to 
trap the second field before it has a chance to roll away, leading to a second epoch of thermal inflation driven by the 
potential 

V = V2-ml\(l)2\'' + ... (41) 

The residual flatons left after parametric resonance from the first epoch of thermal inflation may be troublesome if 
they do not decay before nucleosynthesis. Their abundance evaluated after the second epoch of thermal inflation is 



^ ^ g4TpRiy/*V^^^A2 e2m^,V2 



(42) 



..-14.5 f Ml Y AO^^GeV y feA / Tp, \ 
' lOWGeVy ^ Ah J \e2j \GeY J 



m^J VsOOGeV/ I m2 M f ) [ V2 



2 



(43) 



Thus, a second epoch of thermal inflation may significantly dilute the residual flatons from a first epoch, which could 
remove the restriction Mi < 10^"* GeV which is otherwise demanded by nucleosynthesis. Conceivably one may in 
this way make thermal inflation viable with a GUT Higgs field or even with a modulus, though more investigation is 
needed to see whether this is a real possibility. 

Henceforth we will assume that Mi is sufficiently small to allow (pi to decay before nucleosynthesis, which allows 
us to take M2 small enough to have a comfortably high final reheat temperature. 



F. Solving the moduli problem with double thermal inflation 

For simplicity we take the abundance of big bang moduli before thermal inflation to be n$/s ~ O.l^QMp^^'^m^^^^ 
and assume that parametric resonance leads to effectively complete reheating so that e ~ 0. These assumptions lead 
to the maximum possible moduli abundance. Each epoch of thermal inflation then increases the entropy by a factor 
Ai ApRi ~ V^^'^ /25Tq^. Therefore the abundance of big bang moduli after double thermal inflation is 



'a- 



^ 10-17.5 AO" GeV y / lOi"GeV y /TciV (Tc2 
^ \ Ml J \ M2 J J \ m* 



(44) 



V300GeV/ \ V \ V-? ^ ' 



MpiJ VaoOGeV/ V ^1 / V ^2 

The abundance of thermal inflation moduli produced at the end of the first epoch of thermal inflation evaluated 
after the second epoch of thermal inflation is 
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(46) 



10 



^16.5/ \* /10i°GeVy / T2 



1013 Gcvy ^ M2 J 

\MpJ V300GevH"^|M2y' ^ ^2 / 



The abundance of thermal inflation moduh produced at the end of the second epoch of thermal inflation is 



(48) 



10i"GeVy \MpiJ V "^'i- / \mlM^ 

We see that two independent bouts of thermal inflation can solve the moduli problem for a wide range of the vevs, 
even if parametric resonance is extremely efficient. 



G. Topological defects 

We end this paper with a brief discussion of the cosmological production of topological defects, namely walls, strings, 
monopoles and textures. 

Each type of topological defect is associated with a scalar field (in general multi-component) with nonzero vev. 
Among several possibilities, we consider here only two cases. The first is that the vev belongs to a GUT higgs 
potential, and that it has the non-flat form usually considered. The second is that the vev belongs to a flat potential. 

For a GUT higgs with the standard non-flat potential the temperature after inflation is never high enough for the 



defects to form by the usual Kibble mechanism 33|. (We are not of course concerned with any defects forming 
before ordinary inflation since they have been diluted away.) They can only form near or at the end of ordinary 
inflation, and even that requires that the bound Eq. (^) on the inflationary potential is saturated |26| 



Consider first monopoles, using the standard results |3^. The abundance of monopoles, after some initial annihi- 
lation, settles down soon after the GUT transition to a value n/s ~ 10"!". The strongest bound on their present 
abundance comes from baryon decay catalysis in neutron stars, which requires n/s ^ lO"^^, Thus the entropy must 
increase by a factor 10^^ between the end of ordinary inflation and the present. If reheating after ordinary inflation 
is prompt, the factor is the one A arising from thermal inflation. We see from Eq. ( ^3|) that a single bout of thermal 
inflation is probably insufficient, but two bouts could be enough. Alternatively, if reheating after ordinary inflation is 
long delayed this gives an additional increase Aord ~ 10^^ GcV/Tr, which could be enough to make just one bout of 
thermal inflation viable. 

Depending on the GUT symmetry, gauge strings might also be produced, which would be cosmologically significant 
perhaps providing candidates for the origin of large scale structure. On cosmological scales their evolution is not 
affected by thermal inflation because their spacing is outside the horizon during that epoch. (This is just the statement 
that there are much less than 50 e-folds of thermal inflation.) The same applies to other defect networks formed before 
thermal inflation (global domain walls, monopoles, strings or textures). 

Consider now defects associated with a fiat potential. They form if at all at the end of thermal inflation. Consider 
first the case of Z„ symmetry (Sections 2.1 and 2.3). A discrete symmetry used to be regarded as problematical 
for cosmology, because when it is spontaneously broken it seems to lead to cosmologically forbidden domain walls. 
However, if the symmetry is also explicitly broken, as will typically be the case for the flaton potential, there need be 
no problem because walls do not necessarily form and if they do form they do not necessarily survive (because the 
vacua on either side of a wall may have different energy density). If, on the other hand, it is exact it will probably be 
a discrete gauge symmetry which again avoids the domain wall problem because there is only one physical vacuum. 

If there is a global U{1) symmetry, strings can form at the end of thermal inflation with the strings later joined 
by walls if the symmetry is approximate. An example of this might be Peccei-Quinn symmetry p^ . Local strings 
forming at the end of thermal inflation would have too little energy to be cosmologically signiflcant. Finally, if the 
flaton fleld giving rise to thermal inflation has two or more components as in Section 3.5 then monopoles or textures 
might form at the end of thermal inflation but we have not considered this case. 
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IV. SUMMARY AND CONCLUSION 



Flatons are scalar fields with masses m of order 10^ to 10"^ GeV and vacuum expectation values M ^ m. They 
arise naturally in supersymmetric theories and indeed it is not unreasonable to suppose that they are the only source 
of vevs in this range (in the observable sector). Flatons with M > 10"'^'* GeV are cosmologically dangerous, and in 
particular moduli with M ~ Mpi are overproduced by twenty orders of magnitude in the standard cosmology, which is 
the well-known Polonyi/moduli problem. In this paper we have explained how the problem may be solved by flatons 
with smaller vevs, in the range 10^ GeV ^ M < lO^'^ GeV that is theoretically very natural for flatons other than 
moduli. 

Such flatons solve the moduli problem by generating an era of thermal inflation. Thermal inflation occurs when 
the flaton is held at zero by thermal effects, and it typically lasts for about 10 e-folds and occurs at a very low energy 
scale. These properties are precisely what is required to sufficiently dilute the moduli produced before thermal inflation 
without affecting the density perturbation produced during ordinary inflation (10 e-folds), while not regenerating them 
again afterwards (low energy scale). Detailed calculations show that a single epoch of thermal inflation driven by a 
flaton whose vev is within one or two orders of magnitude of 10^^ GeV can solve the moduli problem, though the 
constraints are quite tight. 

It is easier for thermal inflation to rescue flatons with vev 's M > 10^^ GeV (in particular, moduh with M A/pi) 
if the latter do not themselves give rise to thermal inflation. Remarkably, segregation of flatons into a class which 
thermally inflate and have small vevs, and a class which do not and have large vevs is exactly what one expects from 
a theoretical viewpoint. The larger the vev of a flaton, the less likely it is to be trapped at the origin in the early 
Universe, because the finite temperature contribution to the effective potential becomes ineffective at field values 
bigger than the temperature. 

There are several aspects of cosmology which we have not addressed in the present paper, notably axion cosmology 
and baryogenesis which will be the subjects of respectively [ p^Jl3| . Let us close by briefly discussing the latter topic. 
As successful thermal inflation sufficiently dilutes any pre-existin g m oduli it will also dilute any pre-existing baryon 
number to negligible amounts. However, as will be discussed in ||l^, there are several possibilities for baryogenesis 
within the context of thermal infiation itself. One especially promising mechanism can occur if the flaton which 
gives rise to thermal inflation also generates the mass of a right-handed neutrino. A lepton asymmetry can then 
be generated after thermal inflation. The partial reheat temperature after thermal inflation can be high enough 
to restore the electroweak symmetry, and so this lepton asymmetry can be converted into a baryon asymmetry by 
non-perturbative electroweak effects O]. 



APPENDIX 



To arrive at a more sophisticated estimate of the moduli abundance we solve the equation of motion of the modulus 
in the potential Eq. (E4), which is 



(5$ + 3i?(5$ 



m 



(5$ + aH^ (5$ - $0) = 



(50) 



where we take H = p/t with p — 1/2 for radiation domination and p = 2/3 for matter domination. One would expect 
a ^ 1. An estimate of $0 can be obtained by taking the distance between the self-dual points of the target space 



modular symmetry SL(2,Z) ||5^ using the usual orbifold Kahler metric for the moduli. This gives $q ~ 0.1. In this 
case one can easily check that our approximation of neglecting the contribution of the moduli to the energy density of 
the Universe before the asymptotic solution is attained is consistent. (5$ will rapidly settle to (5$ = $0 when H ^ m$, 
and so we take ^$(0) — $0 and ^<I>(0) — 0. With these initial conditions, Eq. (BQ) has the solution 



(5$ = ap^o 



where s^.j/ is a Lommel function, /i = —3(1 — p)/2 and i'^ 



S<i> 



a^Q (p\ ^ (1 + II 



-ap 



{ip- 1)^/4. At late times 

3p 



H 



sin TO$t 



(2 - 3p)7r 



(51) 



(52) 



The coefficient has a weak dependence on a and p, and for a ~ 1 and p = 1/2 or 2/3 we get to a good approximation 
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^^-3*°Uj + (53) 



Therefore the moduh abundance is given by 

n$ <I>5 



(54) 



Setting p = 1/2 gives Eq. (||) and setting p = 2/3 gives Eq. (p5|). 
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